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ABSTRACT 

We use data from the SDSS to investigate the evolution of the large-scale galaxy bias as a 
function of luminosity for red galaxies. We carefully consider correlation functions of galax- 
ies selected from both photometric and spectroscopic data, and cross-correlations between 
them, to obtain multiple measurements of the large-scale bias. We find, for our most robust 
analyses, a strong increase in bias with luminosity for the most luminous galaxies, an inter- 
mediate regime where bias does not evolve strongly over a range of two magnitudes in galaxy 
luminosity, and no evidence for an upturn in bias for fainter red galaxies. Previous work has 
found an increase in bias to low luminosities that has been widely interpreted as being caused 
by a strong preference for red dwarf galaxies to be satellites in the most massive halos. We 
can recover such an upturn in bias to faint luminosities if we push our measurements to small 
scales, and include galaxy clustering measurements along the line-of-sight, where we expect 
non-linear effects to be the strongest. The results that we expect to be most robust suggest that 
the low luminosity population of red galaxies is not dominated by satellite galaxies occupying 
the most massive haloes. 
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1 INTRODUCTION 

The evolution of dark matter can be accurately modelled within the 
A Cold Dark Matter (ACDM) scenario, giving the number density 
of collapsed objects (e.g., dark matter halos) of given mass as a 
function of redshift. In the widely accepted model of galaxy evolu- 
tion, galaxies form within the gravitational potential wells of host 
dark matter haloes (White & Rees 1978}. Thus, we can reduce the 
problem of modelling galaxy populations to one of working out the 
set of assumptions on how to relate luminous with dark matter. In 
fact, simple instructions can provide an excellent description of the 
observed clustering of galaxies in the near (see, e.g., |Norberg et al.] 
2002, Zehavi et al. 2010) and in the more distant universe (see, e.g., 
|coil et al.|2006||McCracken et al.|2007> . 

In the galaxy formation framework, structure builds hierarchi- 
cally, with smaller halos being the first to collapse and merge to- 
gether to build up larger haloes over time. Galaxies, as permanent 
residents of these haloes, must therefore also grow through merg- 
ing, but the link between galaxy growth (in terms of stellar mass) 
and halo growth (in terms of dark matter mass) is far from direct. 
Firstly, stellar mass in galaxies may also grow through the process 
of turning cold gas into stars, and secondly galaxy growth through 
merging is slower than halo growth - galaxies do not always merge 
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when their host haloes do, and there is strong evidence that stel- 
lar mass may not be conserved during such a process (White et al. 
[2007l|Conroyetal.|2007^ . 

There is a firmly-established colour bimodality in the galaxy 
population seen at low redshift (e.g. |Blanton et al" 2003| >, and this 

This striking 



continues to earlier epochs (e.g, |Bell et al.||2004 



separation of the galaxy population into a blue (star-forming) and 
red (quiescent) clouds must therefore be an outcome of any galaxy 
evolution model. Understanding how these two populations evolve, 
and how galaxies go from one cloud to the other (i.e. what process 
quenches star formation), has been a long-standing quest in galaxy 
evolution, and one for which clustering studies have the potential 
of being particularly insightful. 

The large-scale clustering strength of a particular population 
of galaxies is directly linked to the host dark matter halo mass 
(Bard een et al.|1986|[Cole & Kaiser|1 989). Consequently, study- 
ing the clustering properties of a sample of galaxies, as well as its 
evolution with redshift, remains as one of the most powerful ways 
to constrain galaxy evolution. Many previous studies have inves- 
tigated the clustering of galaxies as a function of their colour and 
luminosity (recent studies include, e.g., Willmer et al. 1998] [Nor^ 
berg et al.|2002| |Madgwick et al.|2003| |Zehavi et al.|2005||Croton 



et al.|2 006 Ros s et al.|2007| |McCracken et al.|2008| |Ross et~ 
2009l|Tinker & Wetzel|2010||Zehavi et al.|2010| >. In broad terms, it 
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has become clear that the large-scale clustering strength increases 
with luminosity, and with redder colour. 

Clustering measurements of galaxy sub-populations, in one 
sense, represent statistical restatements and refinements of the well 
known morphology-density relationship (Dressier 19801, which is 
now known more accurately as a colour-density relationship (see, 
e.g., |Ball et al.|2008[|Skibba et al.|2009") . We now know that the bi- 
modality in the colours of galaxies extends to their large-scale clus- 
tering strengths, revealing that a halo-mass dependent effect is im- 
portant in delineating blue from red galaxies. Essentially, the most 
luminous, most red galaxies occupy the highest peaks in the den- 
sity fields, and for any particular luminosity, red galaxies are found 
to occupy higher peaks than their blue counterparts. This suggests 
red galaxies were the first to begin their hierarchical mass assembly 
— qualitatively consistent with the interpretation that they are red 
because they have used up/been stripped of their cold gas (see, e.g., 
|Larson et al.[T980l|Cowie et al.|1996) . 

If we now focus on the clustering of red galaxies, there is inter- 
esting evidence that their large-scale clustering strength increases 
towards the faint end (see, e.g., Norberg et al.|20 02; Zeha vTet al.| 



|2005HSwanson et al.|2008||Cresswell & Percival|2009) . The natural 
interpretation is that faint red galaxies are predominantly satellite 
galaxies in large mass haloes. This would further imply that red 
galaxies with luminosities ~L* prefer less massive halos than their 
brighter and fainter counterparts — qualitatively consistent with the 
picture of SO galaxies occupying groups and luminous and dwarf 
red galaxies occupying clusters. 

While it is uncontroversial that a significant percentage of faint 
red galaxies are satellites in large mass haloes, the exact proportion 
that are satellites, and whether this is sufficient that these galaxies 
dominate the population, is less clear. Wang et al. ( 2009 ), selecting 
on colour, find that dwarf red galaxies (fainter than those studied 
in any of the other cited works here) are predominantly central, 
rather than satellite galaxies. This is in contrast to Hain es et al.| 
(20071, who, when selecting on Ha emission, find nearly all red 
dwarf galaxies are satellites. Furthermore, there is some tension 
between analyses of large-scale bias and the work of Masje dTet al.| 
(20081, as they find that faint red galaxies have a very weak con- 
tribution to this growth of luminous red galaxies (LRGs), i.e. if the 
faint galaxies occupy the same halos as LRGs, they are somehow 
precluded from merging with them. Finally, Brown et al. ( 2008 ) do 
not see an upturn in the large-scale bias towards fainter luminosities 
when studying red galaxies in the Bootes field, using photometric 
redshifts. 

Motivated partly by the tension described in the previous para- 
graph, and partly by the desire to study fainter red galaxies, in this 
work we take a new look at the large-scale bias of red galaxies in 
the SDSS seventh and final data release, using both photometric 
and spectroscopic data. The photometric data allows us to analyse 
a significantly larger volume of the Universe containing the lower 
luminosity galaxies than the spectroscopic data. This means we can 
consider the same volume to analyse a wide range of galaxy lumi- 
nosities. We also use the cross-correlation between these data with 
brighter spectroscopic data to confirm the robustness of our results. 

This paper is organised as follows: In Section [2] we describe 
the photometric and spectroscopic redshift catalogs we use for our 
measurements. In Section [3] we describe how we measure corre- 
lation functions (angular auto- and cross-correlation functions and 
redshift space 3D auto-correlation functions) and how we use these 
measurements to measure the bias of galaxy samples. In Section]?] 
we present the bias measurements from our photometric redshift 
samples, and in Section[5]from our spectroscopic redshift samples. 



In Section|6] we compare our spectroscopic and photometric results 
to each other and to the results of previous studies. In Section^] we 
discuss the physical implications of our measurements and in Sec- 
tion]^] we present a summary of our conclusions. 

Throughout, we use M r as shorthand for M r — 51og(/i). We 
assume a flat cosmology with Q. m = 0.25, h = 0.7, as = 0.8, and 
fit = 0.045. 



2 DATA 

We use data from the SDSS seventh data release (DR7). This sur- 
vey obtained wide-field CCD photometry i Gunn et al.|19 98) in five 
passbands (u, g, r, i, z; e.g., |Fukugita et al.|1996} , amassing nearly 
10,000 square degrees of imaging data for which object detection is 
reliable to r ~ 22 ( Abazajian et al. 2009). From these photometric 
data, spectroscopic targets have been identified and observed yield- 
ing a sample of galaxies with over 600,000 spectroscopic redshifts 
complete to a Petrosian ( 1976J magnitude limit of r < 17.77 occu- 
pying over 8,000 square degrees (Abazajian et al. 2009). We utilize 
data from both the spectroscopic redshift catalog (as this allows 
precise redshift information) and a photometric redshift catalog (as 
this allows us to go more than 2 magnitudes deeper). 

The photometric redshifts used were obtained from the SDSS 
photoz table. Redshifts in this catalog were estimated using a hy- 
brid template/empirical approach, and the output includes rest- 
frame absolute magnitudes, k-corrections, and galaxy-type values. 
We include all galaxies with z p hot < 0.1 and M r < —17.75 
(which makes our sample approximately volume limited for r < 
20). In order to select red galaxies, we use the (k-corrected) double 
colour cut 



u-r > 2.2, g-r> 0.8, 



(1) 



which ensures that the galaxies are truly red and there are few blue 
interlopers in our sample. Fig.[T]displays the distribution of u — r 
colours of the galaxies in our red sample, with a solid line. The 
dashed displays the distribution we would have had if we simply 
selected u — r > 2.2. Our double cut clearly removes galaxies that 
are close to our cut limit, giving us a distinct sample of galaxies 
(more so than selecting on type value provided in the photoz table). 
We find that performing this double-cut also improves the reliabil- 
ity of the photometric redshift estimates. We further cut the data to 
the angular mask of Ross et al. (2010; hereafter R10), leaving just 
over 6000 square degrees of observing area in the Northern con- 
tiguous area of SDSS, and subdivide this sample into three subsam- 
ples with -17.75 < M r < -18.75 (27023 galaxies), -18.75 < 
M r < -19.75 (23081 galaxies), and -19.75 < M r < -20.75 
(17156 galaxies). 

We also consider galaxies selected from the main SDSS spec- 
troscopic sample. The red galaxy samples were selected using the 
colour cut as defined by Eq. {T}, dividing into five subsamples, 
by absolute magnitude, —17 < M r < —18 (2082 galaxies), 
-18 < M r < -19 (11663 galaxies), -19 < M r < -20 
(54908 galaxies), -20 < M r < -21 (118853 galaxies), and 
-21 < M r < -22 (97947 galaxies). Note that we place no red- 
shift limits on the red galaxy samples we use. 

We also consider spectroscopic samples without colour se- 
lection to cross-correlate against the photometric redshift samples. 
These samples were selected to be approximately volume limited; 
one is limited to z < 0.15, M r < -20.75 (129435 galaxies) and a 
second limited to z < 0.1, M r < -19.75 (140341 galaxies). 
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Figure 1. The distribution of u — r colours of galaxies with z v j %ot < 0.1 
and M r < —17.75 (dotted line) and for galaxies with only u — r > 2.2 
and g — r > 0.8 (solid line). 



3 METHODOLOGY 

We now outline the methodology used to calculate the bias of red 
galaxies as a function of their absolute magnitude. For the photo- 
metric data and for cross-correlations with the spectroscopic data, 
we use the projected angular correlation function. For the spectro- 
scopic data we consider the full 3D correlation function. 



3.1 Calculating Angular Correlation Functions 

In order to calculate angular auto-correlation and cross-correlation 

] estimator 



functions, w(6), we use the Landy & Szalay 1 1993 
D 1 D 2 {9) 



■m = 



DiR 2 (9) - D 2 Ri{9) + RiR2{6) 
R 2 R 2 {9) 



(2) 



where Di and D 2 (and Ri and R 2 ) represent separate data samples 
for the cross-correlation and are the same for the auto-correlation. 
Our random catalog representing the photometric data contains 10 
million points and is the same as used in R10. The photometric cat- 
alogs (data and random) are masked for seeing, reddening, bright 
stars, satellite trails, etc. in the same manner as R10. At scales 
larger than 1°, we employ a pixel based routine using SDSSpix. 
We use the pixel resolution (64) such that the pixels have an area 
with equivalent circular radius of 0.083°. We confirm that the point 
to point and pixelized methods agree at overlapping scales. The 
random catalog representing the spectroscopic data contains over 6 
million points and was created as described in (Reid et al. 
with a pixel based mask based on HE ALPIX (|G6rski et al 

We use a jackknife method (e.g., |Scranton et al.|2002) , with 
inverse-variance weighting to estimate the errors and covariance 
matrix (e.g., Myers et al.||2007 i for w(9). The method is nearly 
identical to the method described in detail in Ross et al. (2007; here- 
after R07) and applied to DR7 data in R10. The jackknife method 
works by creating many subsamples of the entire data set, each with 
a small part of the total area removed. R07 found that 20 jack-knife 
subsamplings are sufficient to create a stable covariance matrix, and 
we therefore use 20. These 20 subsamples are created by extract- 



ing a contiguous grouping of l/20th of the unmasked pixels in 20 
separate areas. Our covariance matrix, C w , is thus given by 



(3) 



= H Z)fclit w /«»( i) - Uk(0i)][ufuu(0j) - ^k{9j)}, 

where uik{9) is the value for the correlation measurement omitting 
the fcth subsample of data and i and j refer to the i th and j angu- 
lar bin. The jackknife errors are simply the square-root of diagonal 
elements of the covariance matrix. 



3.2 Calculating 3D Correlation Functions 

In order to calculate the 3D auto-correlation function, £(r, /i), we 
again employ the Landy & Szalay ( 1993 ) estimator but in this case 
DD, DR, and RR are now functions of r, the physical separation 
we calculate given our assumed cosmology, and also /i, the cosine 
of the angle between the radial direction and the alignment of a 
galaxy pair. Further, we must assign the randoms radial positions, 
which we do by sampling a fine resolution cubic spline fit to the 
galaxy redshift distribution. We also weight each galaxy according 
to the number density at its radial position (Feld man et al.|1994] >: 
this weight is optimal if the galaxies Poisson sample the density 
field. 

We analytically estimate the errors and covariance of the 
£(r, n) measurements. Given optimal weights, w(r) = 1/(1 + 
n(r)P(k)), and the power spectrum, P(k), its error can be esti- 
mated by (Tegmark 1997): 

(27T) 3 P(fc) 2 



V k V eff (k) 



where V e ff(k) = J 



n(r)P(k) 



d s r and V k 



(4) 



4ivk 2 Ak. 



l+n(r)P(k) 

Given that the correlation function is just the Fourier transform 
of the power-spectrum, we can estimate the covariance, C(, in the 
measured £ between bins centred at ri and r 2 as 



Q(ri,ra) = 



1 

2^2 



dkP(k) 2 sm(rik)sin(r 2 k) 
r ir2 V eff (k) 



+ Sr 1 ,r 2 /n p , (5) 



where 8 is the Rronecker delta, and n v is the number of pairs in- 
cluded in the bin. Note that although this procedure ignores corre- 
lations induced by the sample window function, it should be suffi- 
ciently accurate for the analysis attempted in this work, as we only 
wish to measure the broad amplitude of the clustering rather than 
the relative clustering on different scales: for such measurements 
of amplitude, correlations between data are unimportant. 



3.3 Theoretical Modelling 

We model the non-linear power-spectrum using the fitting formulae 
of |Smith et allpOOl} , which we use in combination with the trans- 
fer function of Eisen stein & Hu| < | 1998) to include the effects of 
baryons. Given this power-spectrum, we determine the isotropic 3- 
dimensional real-space correlation function £(r) via Fourier trans- 
form. We model the redshift-space correlation function as ( |Hamil-| 
|ton|1992| > 



where 

&(r) 
Mr) 



(b 2 + lbf+lfnr), 
(f&/+y/ 2 )K(r)-£'(r)], 



(6) 

(7) 
(8) 
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ii('"> = ^/ 2 K(r) + ^'(r)-^"W], (9) 



In the regime where linear theory is a good approximation this re- 



Pt are the standard Legendre polynomials, and 



^" = 5r- 5 / e(r')(r') 4 dr', 



(10) 



(11) 



b is the large-scale bias of the galaxy population being considered. 
We will be comparing our models to measurements limited to dif- 
ferent maximum fj, values. Thus, the model will become 



(12) 



In order to calculate model w(8), we must project £ s (/it, r) 
over the radial distribution of galaxy pairs in a particular sample 
(or samples in the case of cross-correlations). 

w(9) — J dzi J dz 2 n i (z 1 )n j (z2)£ s [fi, r(9, z\, z 2 )] , (13) 

where m is the normalized redshift distribution of sample i (and 
i = j for the auto-correlation). The galaxy separation r is a func- 
tion of the angular separation of the galaxies 9 and their redshifts 
z\ and Z2 (as is /i). 

3.4 Modelling Redshift Distributions 

The photometric redshifts we use have significant error associated 
with each redshift estimate. This implies that estimation of the true 
redshift distribution of the photometric redshift data is not trivial. In 
order to estimate the true redshift distribution for each of our pho- 
tometric samples, we follow the procedure outlined in RIO. Given 
a redshift estimate and its error, we assume a Gaussian PDF. We 
then convolve this PDF with the luminosity function and volume 
element/redshift relationship to obtain a final PDF, which we then 
sample 10 times. In this work, we find that adding a bias to the pho- 
tometric redshifts (simply a shift in the mean) of +0.005 for each 
sample allows for the best agreement between all of our clustering 
measurements (see Section[4]l. 

We can construct M r distributions in a similar manner. Given 
the M r of the redshift estimate, we can calculate the value of M T 
at any point along the redshift PDF based on the change in the 
distance modulus, thus obtaining a PDF of M r for each galaxy. 
This allows us to determine the expected mean M r for any particu- 
lar photometric redshift sample. When we calculate angular cross- 
correlations between photometric and spectroscopic samples, the 
redshift ranges of the samples are not matched. Thus, the distribu- 
tion of luminosities of the galaxies that contribute to the measured 
clustering signal will not necessarily match that of either sample. 
In order to find the expected mean M r of galaxies contributing to 
the cross-correlation, we can simply ignore the parts of the redshift 
PDF outside the spectroscopic redshift bound and thereby deter- 
mine the effective mean M r of photometric galaxies between any 
chosen spectroscopic redshift bounds. 

3.5 Calculating Bias Values 

We assume a simple linear bias model between the the 2-pt cluster- 
ing in the dark matter and galaxy over-density fields, such that the 
bias b, is defined by 



P(k)g = b 2 P(k) DM - 



(14) 



lationship naturally extends to £(r) and w(6). The use of the Smith 
|et al.|{2003| > non-linear P(k) has the potential to extend this rela 
tionship to weakly non-linear scales. Thus, for arbitrary limits (but 
always assuming b is scale independent in Eq. | I4[ i we determine 
the best-fit value of b for a particular galaxy sample by minimizing 

x = Kfc) - t&(0,)] czHOiM) [w(e 3 ) - w{e 3 )] (15) 



where the model correlation function w and x are implicitly de- 
pendent on the value of b. For the 3D correlation functions, w(9) is 
replaced by f s (r, /J, max ) and 9 by r. 



4 BIAS FROM PHOTOMETRIC DATA 

We measure the bias of galaxies in our photometric samples by 
measuring their respective auto-correlations and cross-correlations 
with the spectroscopic sample of galaxies with z spBC < 0.15 and 
M r < —20.75 (which we denote "full"). The redshift limit means 
this sample covers most of the redshift distribution of the photo- 
metric samples and the absolute magnitude limit ensures that the 
bias of the galaxies in the sample should not change significantly 
with redshift. Based on w(9) measurements with 9 > 1°, the bias 
of the full spectroscopic sample is 1.32 ± 0.04, which agrees with 
the results of Zehavi et al. ( 20 1 1, as it lies between the bias values 
of their M r < -20.5 and M r < -21 samples. 

We determine the mean M r of each photometric sample by 
estimating the distribution of M r values as described in section 
|3.4| The cross-correlation with the spectroscopic sample implies a 
hard limit of z < 0.15 for the photometric galaxies contributing 
signal to the measured cross-correlations. We find that imposing 
this limit when calculating the magnitude distribution changes the 
mean M r value by less than 0.05 magnitudes in every case. This 
implies that the z < 0.15 limit makes only a small difference to 
any of the magnitude distributions of the photometric galaxies that 
contribute to the cross-correlation signal. Therefore, it is valid to 
expect the bias of the photometric galaxies contributing to the auto- 
correlation signal to be the same as the bias of the photometric 
galaxies contributing signal to the cross-correlation with the z < 
0.15 spectroscopic sample. 

Fig. [5] displays the measured w(9) for our three main photo- 
metric samples. The auto-correlations are shown by the black tri- 
angles and the cross-correlations are displayed in red-squares, with 
luminosity increasing from bottom to top. The solid curves display 
the best-fit model w(9) fit for 1° < 9 < 5°, where we have jointly 
fit the w(9) auto- and cross-correlations (note, bias values can be 
independently determined for the auto- and cross-correlations; the 
joint-fit determines the most consistent value and reduces the over- 
all uncertainty). We restrict the fits to 9 < 5° because systematic 
effects in the data (due to, e.g., reddening and star/galaxy separa- 
tion) and the modelling (effects of redshift space distortions) be- 
come more significant at larger angular scales. The jack-knife er- 
rors increase only slightly over our range of measurements, due 
to competing effects of cosmic variance and shot-noise. The mean 
M r and best-fit b values corresponding to these joint fits are -18.46, 
1.2±0.05, -19.46, 1.43±0.07, and -20.39, 1.40±0.07, with respec- 
tive x 2 values of 13.3, 16.0, and 17.0 (and 13 degrees of freedom). 
We note that these values remain constant (within the la error-bars) 
regardless of the minimum scale we fit to between 0.5° and 2.0°. 

In order to test the robustness of these results to the photo- 
metric redshift determination, we have also obtained estimates of b 
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Figure 2. The measured angular auto-correlations for the three specified photometric samples, all with 2phot < 0. 1 (black triangles) and their cross-correlations 
(red squares) with the spectroscopic sample with z < 0.15, M r < —21.75. The black and red curves display the respective best-fit models (all fit between 1.0 
and 5.0 degrees). 



for galaxies of different luminosity by cross-correlating the photo- 
metric samples with spectroscopic samples covering different red- 
shift ranges. The cross-correlation signal depends on the number of 
galaxies in the photometric sample that lie within the redshift range 
covered by the spectroscopic sample, so testing whether these bias 
estimates match those from our full sample tests the assumed dis- 
tribution of photometric redshifts. We cross-correlate each photo- 
metric sample with two additional spectroscopic samples; one with 
z < 0.1, M T < —19.75 (which denote "near") and second with 
0.1 < z < 0.15, M r < -20.75 (which we denote "far"). We 



find that these redshift limits increase/decrease the mean value of 
M r of the galaxies contributing to the cross-correlation signal by 
~ 0.4 magnitudes. We find the bias of the near sample is 1.17 ± 
0.05 (again in agreement with Z ehavi et al.|2010| and we continue 
to use b = 1.32 ± 0.04 for the far sample. Note that the near and 
far bias values can be determined only from the cross-correlations, 
since auto-correlations of the photometric data cannot replicate the 
effects of a hard redshift cut. 

We display the best-fit b against M r for galaxies in our photo- 
metric redshift samples in Fig. [3] Values determined from the near 
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Figure 3. The measured bias values using the combination of the photo- 
metric auto-correlation and the cross-correlations with the z spec < 0.15, 
M r < —20.75 sample (solid triangles), the cross-correlation with the 
Zspec < 0.1, M r < —19.75 ("near"; open triangles), and the cross- 
correlation with the 0.1 < z spec < 0.15, M r < —20.75 data ("far"; 
open circles). 



Figure 5. The measured bias of red spectroscopic galaxies as a function of 
M r , for four different choices of r m i n and /imax- The M r values have 
each been shifted slightly for clarity. For comparison, the bias values of the 
red photometric galaxies (same as plotted in Fig.[3) are displayed with cyan 
open circles. 



cross-correlation are shown by the open triangles, while values de- 
termined from the far cross-correlation are shown by the open cir- 
cles. The b values determined from the joint fit to the photometric 
auto-correlations and cross-correlations with the full spectroscopic 
sample are plotted with solid triangles. Note that the joint fit esti- 
mates of b serve as our primary results, while the additional cross- 
correlations serve as a (correlated) consistency check. 

The fact that the results from all cross-correlations agree is im- 
portant confirmation that photometric redshift distributions we have 
estimated are accurate. We note that if we do not include the pho- 
tometric redshift bias of +0.005, the cross-correlations do not yield 
consistent results. For every sample, the cross-correlation with the 
near data would imply a much lower (by ~ 30%) value of b and the 
cross-correlation with the far sample would yield a much higher 
(by ~ 25%) value of b, implying (the obviously incorrect result) 
that b(M r ) would be highly oscillatory. However, this photometric 
redshift bias does not have a strong effect on the best-fit b values 
of the three main samples. When the photometric redshift bias is 
set to zero, the best-fit b values decrease slightly to 1.19 ± 0.05, 
1.39 ± 0.07, and 1.38 ± 0.07 with increasing luminosity. 

We have tested our results further by making changes both to 
the data we use and the modelling. We find that our results do not 
change (within la) regardless of the minimum angular scale we 
use to find the best-fit data. We have investigated multiplying the 
error on the photometric redshifts by a constant factor and leav- 
ing this parameter free when finding the best-fit model. We find 
that this can slightly improve the \ 2 values, but not to the extent 
demanded by removing a degree of freedom. Importantly, doing 
so does not change the trend we find with M r — the bias of the 
lowest luminosity red galaxies is always lowest. We find the same 
behaviour with different colour cuts and also if we make changes 
to the redshift/magnitude limits. No matter the analysis we apply, 
we find the same result — red galaxies with M r > — 19 have lower 
bias than their brighter counterparts. The bias measurements of the 



photometric galaxies thus suggest that the bias of red galaxies is 
approximately constant between —19 > M r > 21 and decreases 
to lower luminosities. 



5 BIAS FROM SPECTROSCOPIC DATA 

We measure £, s (fi < fimax, r) in order to determine best-fit val- 
ues of b for the five different spectroscopic samples of red galax- 
ies with magnitudes between —17 and —22. The results of these 
measurements, for Umax = 1 (bottom) and for Umax = 0.5 
(top) are displayed in Fig. [4] The two most luminous samples 
(-21 < M r < -22, magenta; -20 < M r < -21, green) 
display smooth shapes over the range of the plot and are consis- 
tent with respect to each other for both fi ma x values. The least- 
luminous bin (—17 < M r < —18) is quite noisy and displays a 
dramatic difference when Umax is changed — its clustering am- 
plitudes are roughly double those of the most luminous sample for 
r < 10/i _1 Mpc when fimax = 1 but this behaviour is removed 
for Umax = 0.5. The two samples — 18 < M r — 19 (red) and 
— 19 < M r < —20 also display significant dependence on the 
scale and Umax, as both display (relatively) larger amplitudes at 
smaller scales and greater Umax values. We note, however, that any 
cut that is more restrictive than fi < 0.5 does not produce a statis- 
tically significant change in any of our bias measurements. 

Not surprisingly, we find there is a strong dependence both 
on the minimum r and maximum fi that we use when we esti- 
mate b for the spectroscopic samples. Fig.|5]plots b versus M r for 
four separate r m in/ Umax limits. For all bias estimates, the maxi- 
mum scale we fit to is 40/i _1 Mpc. If we use all angles to the line 
of sight and fit above 6/i~ 1 Mpc (black triangles), we find that the 
bias increases as the luminosity of the red galaxies decreases. This 
behaviour is similar to that reported previously by Swanson et al. 
([2008) and also |Cresswell & PercivaT| p009l >. However, if we fit 
for r m in > 12h~ Mpc (open red triangles) or fi < 0.5 (green 



© 2010 RAS, MNRAS 000,[T|frT| 



Large-scale clustering of faint red galaxies 7 



200 




1 10 

r (h^ 1 Mpc) 



Figure 4. The measured redshift-space correlation function multiplied by r 2 for the rive spectroscopic samples: — 17 < M r < — 18 (black), — 18 < M r — 19 
(red), —19 < M r — 20 (blue), —20 < M r < —21 (green), and —21 < M r < —22 (magenta). The top panel displays the measurements when we restrict 
fi < 0.5 and the bottom panel displays the results with no such restrictions. 



squares), we see only weak evidence of an increase in the bias to- 
wards lower luminosity. If we require both that r m in > 12/i _1 
and fi < 0.5 (open blue squares), we recover a monotonic in- 
crease in the bias with luminosity. The best-fit bias values for 
r > 12/i _1 Mpc and fi 0.5 are consistent with the trend we 
found in the photometric data (plotted in Fig. [5] with open cyan cir- 
cles). 

The bias results that we expect to have the least systematic ef- 
fects are those calculated from measurements at the largest scales 
(since we expect linear theory to be most valid at larger physical 



scales) and those calculated for pairs furthest from the line of sight, 
since they should be the least dependent on redshift space distor- 
tions (RSD). Note that we minimize the effects of linear RSD. We 
do allow for linear RSD in our modelling, but studies (see, e.g., |Per-| 
|cival & Wh ite 2009; Jennings e t al.|20ll)) have shown that linear 
RSD models begin to fail at larger scales than real-space clustering 
models. Perhaps more importantly, we minimize any of Fingers-of- 
God (FoG) effects ( Jackson 1972). The effects of FoG imply that, 
by restricting our analysis to transverse pairs, we are preferentially 
discounting pairs of galaxies within dark-matter halos. Note that 
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Figure 6. The measured bias of red spectroscopic galaxies as a function of 
M r , for r > I0h~ Mpc, (i ma , = 1 and r > 12/i -1 Mpc, fimax = 
0.5. Also plotted are the recent results of three other recent studies of the 
clustering of red galaxies. 



for the purposes of measuring the linear bias, this is entirely ap- 
propriate, as the clustering of galaxies within dark matter halos is 
clearly non-linear. 

We thus believe that, of the data included in Fig. [5] the r > 
12/i _1 Mpc, fi ^ 0.5 results are the most trustworthy. The x 2 /DOF 
values are all much smaller for the r > 12/i -1 Mpc fits, as for 
fi ^ 1, the x 2 /DOF values are all greater (and as high as 6) than 2.4 
for r > 6ft -1 Mpc and they are all less than 1 for r > 12ft -1 Mpc. 
For fi ^ 0.5, the difference is not as extreme, but in every case 
the x 2 ^DOF values are reduced when the fit is performed at larger 
scales. For the faintest bin, the bias we measure is sensitive to the 
particular fimax, r min we choose (it continues to vary significantly 
if we increase r m in beyond 12h~ Mpc), suggesting that the sys- 
tematic error associated with this bias measurement is much larger 
than the statistical error displayed in Fig.[5]and subsequent Figures. 



6 COMPARISON TO PREVIOUS RESULTS 

Our most robust results suggest there is monotonic increase in the 
bias of red galaxies with luminosity between —17.5 > M r > 
— 19.5. This disagrees with recent findings. A comparison between 
a selection of our results for the spectroscopic data compared to 
some previously published results is displayed in Fig. [6] We have 
normalized the bias values of previous studies such that the bias of 
red galaxies with magnitudes closest to M r = —20.5 is equal to 
1.4. For the Zehavi et al. 2010| results, we have used the w p (r p ) 
measurements they list in Table B9 at scales 6 < r p < 17/i _1 



Mpc. We simply determine the inverse-variance weighted mean of 
the ratio between each of their samples an the — 20 > M r > — 21 
sample and then take its square root to estimate the relative bias of 
each of their samples. 

All of the results plotted in Fig.|6]agree quite well for the most 
luminous red galaxies with M r < —20. In fact, for M r < —19, 
the agreement between our spectrosco pic analysis with r m i n = 
12/i _1 Mpc, Umax = 0.5 (blue squares), Swanson et al. (2008 i, and 



|Zehavi et al.| ( |2010) is excellent. For less luminous red galaxies, the 
bias is strongly dependent on the sample and method used. We find 
that a choice of r min = 10/i' 1 Mpc, Umax = 1 (black triangles) 
for the spectroscopic data analysis produces a result that is quite 
similar to Cresswell & Percival| 1 12009) and |Swanson et al.| (2008), 
while the fit with r m in = 12h~ Mpc, Umax = 0.5 (blue squares) 
shows strong disagreement. 

The dramatic change in the best-fit bias as a function of r m i n 
and Umax suggests that non-linear effects have a great influence on 
the clustering of faint red galaxies, even to scales ~ 10/i -1 Mpc. 
The correlation function measurements we perform on the spectro- 
scopic samples (see Fig. [4} are consistent with the results of |Hogg| 
|et al.H2003) and |Blanton et al.| ( |2005) ), which found bright and faint 
red galaxies to have the highest over-densities on 8/i~ Mpc scales. 
|Cresswell & Percival| ( |2009") fit their power spectra measurements 
to similarly small scales, as they used k < 0.21/tMpc~ . Fur- 
ther, the measurements of Norberg et al. (2002) were restricted to 
r < 5ft~ 1 Mpc, and their bias estimates are therefore quite differ- 
ent than ours. Larger scales appear to be required in order to escape 
the influence of non-linearities. Swan son et al.| j2008) use a larger 
physical scale (20ft -1 Mpc) for the radius of their cells, but being 
that they use counts-in-cells, the counts in large cells are highly 
covariant with those of smaller scales, and non-linearities may per- 
sist. Further, while all of the studies we cite do attempt to account 
for redshift-space distortions, none explicitly remove pairs oriented 
along the line-of-sight, suggesting the possibility that systematic 
effects may still be present. 

If we focus only on DR7 data with minimal contribution from 
clustering along the line of sight, the result is shown in Fig. [7] (the 
w p (r p ) and w(8) measurements should be dominated by trans- 
verse pairs). Within the la error-bars, only the measured bias of 
the M r ~ —19.5 photometric sample appears inconsistent with 
the rest of the data, but it is consistent to within 2a. Given the 
number of data points (19), this is reasonable statistically. The con- 
sistency of the data can be further confirmed by the fact that a fit 
to b{M r ) = (a(M r + 20)) 3 + b a has minimum \ =7.2 for 
a = —0.35 and b — 1.36 (this suggests only that this data is 
consistent with a smoothly varying function and do not believe it 
should be extrapolated outside of the M r values we present). The 
combination of data suggests that the bias of red galaxies increases 
dramatically for galaxies more luminous than -21.5, the bias of 
red galaxies is nearly constant for —19 > M r > —21, but for 
M r > — 19 the bias decreases significantly. 

A reasonable interpretation of the compiled results is that the 
previously reported increase in the bias of red galaxies as the lumi- 
nosity decreases is correct only for non-linear scales and therefore 
not directly related to the average mass of the halos that faint red 
galaxies occupy. Non-linear effects were made stronger due to red- 
shift distortions and non-linear clustering, further exacerbated by 
the very small volume occupied by the low luminosity galaxies. 
The studies that utilize larger volumes (the photometric data) or 
use results only where the influence of clustering along the line of 
sight is minimized (the spectroscopic data herein with n < 0.5, |Li| 
|et al.|2006[|Zehavi et al.|2010| > agree that the bias of red galaxies 
displays a monotonic increase with luminosity. 



7 PHYSICAL INTERPRETATION 

Our results suggest that redshift distortions and non-linear cluster- 
ing have a large effect on the measured clustering of faint red galax- 
ies, even at large scales. Physically, this makes sense. We know that 
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Figure 7. The measured bias of red galaxies, considering only DR7 mea- 
surements with minimal contribution from pairs along the line of sight. 

a significant percentage of faint red galaxies are satellites in high 
mass halos. This implies that a significant proportion of the pairs 
of faint red galaxies are satellites within the same dark matter halo 
and their true physical separations are as large as a few /i _1 Mpc 
(the Virial radius of a W 14 M Q /h halo is ~ 1 ft _1 Mpc). Thus, even 
given the true 3D distribution of galaxies, we would expect to find 
non-linear effects have more influence on the clustering of faint red 
galaxies than for other galaxy sub-populations. 

Finger-of-God effects exacerbate the the non-linear effects. 
Consider a pair of satellite galaxies oriented perpendicularly to the 
line of sight with a separation of 2ft~ 1 Mpc: Assuming they orbit 
the centre of mass of their dark matter halo at 500km/s, we would 
interpret their physical separation to be greater than lOft^Mpc. 
Thus, in redshift space, FoG effects project non-linear clustering to 
larger scales, and the magnitude of this effect should be largest for 
faint red galaxies. By ignoring pairs close to the line-of-sight, FoG 
effects mean that satellite pairs with apparent large-separations are 
preferentially removed. However, these pairs do not truly have 
large-separations, and should not contribute to the true large-scale 
clustering. This explains why removing these pairs brings the re- 
sults closer in amplitude to the projected clustering results. Given 
that studies such as Li et al. ( 2006 1 have found that the pairwise- 
velocity-dispersion is greatest for red galaxies with M T < —19, the 
FoG effect has its strongest impact on the clustering we measure for 
the faintest galaxies. 

Previous studies have found an upturn in the large scale bias 
of red galaxies towards faint luminosities, and this behaviour has 
been attributed to the fact that these galaxies are predominantly 
satelli tes in larger mass haloes (|Blanton et a l. 2005 ; Swanso net al.| 
[2008] |Cresswell & Percival|2009) . In Section [T] we partially mo- 
tivated this work by pointing out that the implications of Masjedi 
et al. (20081 and Wang et al. (20091 suggested a possible tension 
with the previous clustering studies. The results we present in this 
paper show that the reports of an upturn in the bias relationship 
likely cannot be interpreted as an upturn in the linear bias, due to 
systematic effects. We find that the large-scale, linear, bias of red 
galaxies monotonically decreases with decreasing luminosity. This 
in turn suggests that satellite faint red galaxies, residing in massive 



haloes, are not dominant by number with respect to central faint red 
galaxies or faint red satellites of ~ L* galaxies. 

Measurements of the satellite fraction in the red cloud vary. 
Tinker & Wetzel ( 20 1 1 suggest a value of 30% which does not vary 
significantly between z = 0.4 and z — 2. Brown et al. (2008| find 
that the fraction of stellar mass in satellites increases with host dark 
matter halo mass, and such that in haloes with Mu > 1O 14 M0 jh 
there is more mass in satellites than centrals. At z=0, they conclude 
that 30% of the stellar mass in the red sequence (not number of 
galaxies) is in satellites. Wang et al. ( 2009 1, using a group catalogue 
and studying very faint (dwarf) red galaxies, conclude that at only 
around 45% are satellites. 

|Zehavi et al.| ( |2010") suggest an incredibly large satellite frac- 
tion of 90% at the faintest end (-19 < M r < -20), up from 
33% at the intermediate luminosity range (—21 < M r < —20). 
Interestingly, the HOD model that gives these values (driven by the 
small-scale clustering) over-predicts the large-scale clustering am- 
plitude. An alternative HOD model, that only allows red satellites 
in haloes with Mh > 10 13 Mq / h, reduces the total satellite frac- 
tion in the red cloud to 34% and provides a significantly better fit 
to the data at large scales. Note the apparent contradiction to previ- 
ous studies; forcing satellite red galaxies to only occupy high mass 
halos reduces the large scale bias, due to the fact that such a model 
allows a lower satellite fraction. Combined, the results of previous 
studies suggest that the satellite fraction of red cloud galaxies is 
35% and may increase towards the faint end. This suggests that, 
while the satellite fraction allows the bias of red galaxies to remain 
high relative to blue galaxies, the propensity of faint red galaxies to 
be central galaxies in low mass halos allows their bias to decrease 
monotonically. 

Central faint red galaxies present an interesting puzzle in 
galaxy formation. It has been observationally established that the 
satellite population is in general redder that central galaxies of the 
same mass (see, e.g., |van den Bosch et al.|2 008 Yang et al .|2008[ 
|Guo et al.|20l0l >. This can be explained by environmental effects 
as the galaxy falls onto a more massive halo: a slow or a sudden 
stripping of the gas reservoir, due to infall, ram-pressure or tidal 
heating (see, e.g,.|Gunn & Gott|1972[ [Larson et al.|1980[|Moore et| 
al. 1996 ; Balogh & Morris 2000 1. However, a central galaxy has not 
experienced a direct infall onto a larger halo. One possibility is that 
put forward by Wang et al. (20091, in which galaxies may experi- 
ence some form of these environmental effects when passing by a 
dark matter halo (within 3 times the virial radius), without falling 
in. Another possibility is simply that the process that shut down 
star formation in these galaxies in different, or driven by feedback 
effects. 



8 SUMMARY AND CONCLUSIONS 

We have computed the large scale bias of red galaxies in the SDSS 
DR7 as a function of luminosity. We use both photometric and 
spectroscopic data to obtain robust results and utilize the largest 
volume possible. Our major findings are: 

• The bias of red galaxies we measure in our photometric data im- 
plies the faintest red galaxies have the lowest bias (see Fig. [3}. This 
result is independent of any treatment we apply to the data or its 
analysis. 

• The bias increases dramatically with increasing luminosity for 
galaxies with M r < —21. There is no evidence for evolution in the 
bias for galaxies with — 19 > M r > —21, and there is weak evi- 
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dence for a decrease in bias with decreasing luminosity for fainter 
galaxies (see Figs. 5, 6, and 7). 

• The bias of red galaxies measured in our spectroscopic samples 
depends greatly on how we choose to analyze the data (see Figs. 
5 and 6). The treatment that we believe is most free of systematics 
(fmin > 12/i~ 1 Mpc, fi max < 0.5) recovers amonotonic increase 
in the bias of red galaxies as a function of luminosity. This sug- 
gests that previous reports of an upturn at low luminosities in the 
large-scale bias red galaxies were systematically affected by a com- 
bination of redshift-space distortions and non-linearities. 

The results that we show here, based on the large-scale clus- 
tering, give only an incomplete picture of what is happening at the 
faint end as we cannot disentangle the contribution from central 
and satellite galaxies to the overall observed bias value. In some 
sense the bias of faint red galaxies represents the interplay between 
the number of satellites (with strong bias, in high mass haloes), and 
centrals (with low bias, in low mass haloes). Nonetheless our re- 
sults are in good qualitative agreement with other measurements of 
the fraction of red satellites in the sense that the fraction of satellites 
is not expected to be dominant even at the faint end. Treatment of 
the clustering of red galaxies utilizing halo-occupation-distribution 
modelling to constrain the evolution and mass-assembly of red 
galaxies will be presented in a follow-up paper. 
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